Collective oscillations in ultra-cold atomic gas 
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Using both fluid and kinetic descriptions, where repulsive forces between near by atoms are 
included, we discuss the basic oscillations and waves of a cloud of ultra-cold atoms confined in a 
magneto-optical trap. The existence of a hybrid mode, with properties similar to both plasma and 
acoustic waves is described in detail. Tonks-Dattner resonances for confined hybrid modes in a 
spherical cloud are discussed and the prediction of a nonlinear coupling between the dipole resonanc 
and the hybrid modes is considered. Landau damping processes and quasi-linear diffusion in velocity 
space are also discussed. 

PACS numbers: 



I. INTRODUCTION 

In recent years, and mainly due to the emergence of 
laser cooling techniques 1 there has been an increasing 
interest on the physics of ultra-cold atoms. This interest 
has been mainly concentrated on the study of Bose Ein- 
stein condensates, and on the theoretical and experimen- 
tal understanding of their properties [3] . However, at- 
tention has also been given recently to the study of collec- 
tive oscillations of non-condensed atomic molasses in the 
magneto-optical traps |U [51 [B]. From this work, collec- 
tive behavior similar to that observed in plasma physics, 
has emerged, leading to the discovery of an equivalent 
electric charge for the neutral atoms, and of electrostatic 
type of interactions between nearby atoms [7j, and to 
Coulomb like explosions of the atomic cloud when the 
magnetic confinement is switched off [8] . The theoretical 
modeling of such collective processes is still not well es- 
tablished, and this is the main motivation of the present 
work. Here we propose to apply the well known methods 
of waves and oscillations in plasmas to such processes in 
a neutral gas. 

In this work we consider the collective behavior of 
ultra-cold atom gas, in order to identify the basic mech- 
anisms of oscillation, to establish its frequencies and to 
derive dispersion relations for its basic modes of propa- 
gation phenomena. We will use both fluid and kinetic 
descriptions where the main forces associated with the 
laser cooling processes are retained. 

This paper is organized as follows: in Sec. II, we de- 
fine some essential parameters to describe the features of 
a cold atomic cloud confined in a magneto-optical trap. 
The basic set of fluid equations is also established. In 
Sec. Ill, we report the existence of brand new hybrid 
mode for short wavelength scales, somewhere between 
plasma waves and acoustic waves. In Sec. IV, we adress 
the oscillations for the long wavelength scales, where we 
discover modified Tonks-Dattner resonances, which cor- 
respond to confined hybrid oscillations inside the atomic 
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cloud, formally similar to those in ref. [TT]. Nonlinear 
hybrid resonances driven by the dipolar oscillations are 
predicted in Sec. V. We also use a wave kinetic approach, 
in order to refine the fluid description and describe new 
phenomena. Hence, in Sec. VI we describe the system 
via a Fokker-Planck equation that can be directly derived 
from the Schroedinger equation for the collective field of 
the atoms in the radiation and trapping fields. It is shown 
that the waves and oscillations in the gas can be damped 
by a resonant atomic interaction with the collective oscil- 
lations, which is a manifestation of Landau damping. We 
will also show that a spectrum of collective oscillations 
will lead to diffusion in the atomic velocity space, thus 
preventing the laser cooling process to proceed further. 
This effect is added to other diffusion processes already 
known in the literature. Finally, some conclusions are 
stated. 



II. BASIC EQUATIONS AND FLUID 
DESCRIPTION 



The simplest model to describe a gas of cold atoms in a 
magneto-optical trap is based on Doppler cooling jSJ ITU] 
and on the spatial confinement due to the presence of 
a magnetic field gradient. The relevant average forces 
acting on a single atom are based on the quasi-resonant 
radiation pressure force and can be written (for a con- 
venient choice of the relative polarisations of the laser 
beams and the atomic transition) for each direction 
(n = x, r 2 = y, r 3 = z) as 
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This expression relies on the low intensity Doppler model 
for the magneto-optical force (incident on-resonance sat- 
uration parameter per beam Sj„ c = Io/I S at ^ 1)- The 
Zeeman shifts (described by ^r,) and Doppler shifts 
(kLVi) are responsible for trapping and cooling respec- 
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tively. Here Iq is the laser intensity of laser beams inci- 
dent along the six directions, T the natural line width of 
the transition used in the cooling process and A the fre- 
quency detuning between the laser frequency ujl — k^c 
and the atomic transition frequency io a t- Assuming sym- 
metric forces (fj,i = /i) along each of the three directions 
Ox, Oy and Oz one can write to first order in r and v 



atomic clouds and due to the laser cooling beams, cor- 
respond to the simplest possible description of the laser 
cloud interaction, and can be used in a first approxima- 
tion to model the fluid dynamics of the ultra-cold gas, 
which can be derived by computing the zeroth and the 
first momenta of the Fokker-Planck equation, neglecting 
the diffusion term [5]. The basic set of equations can 
then be written as 
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where k is the spring constant of the trap and a the 
friction coefficient, 



a = —8hk^ s 
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which is related to the dipole frequency 
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where M represents the mass of a single atom. Even 
though most magneto-optical traps using for instance al- 
kaline atoms (such as Na, Li, K, Rb or Cs) are not com- 
pletely described by such a simple Doppler model [T] , the 
forces above constitute a good first approach to describe 
the dynamics of single atoms in a magneto-optical trap. 

This description of a magneto-optical trap is known to 
be limited to only moderate number of atoms (typically 
10 5 ). For larger atom number additional forces need to be 
taken into consideration. A second force to be considered 
is the shadow force, or absorption force, Fa, and was 
first discussed by Dalibard [TS]. This is associated with 
the gradient of the incident laser intensity due to laser 
absorption by the atomic cloud. It is an attractive force 
which can be determined by a Poisson type of equation, 
as 



V • F A = - 



n{f) 



(5) 



where n(r) is the atom density and gt the laser absorp- 
tion cross section. Finally, a third force, Fr, can be called 
repulsive force or radiation trapping force, and was first 
considered by Sesko et al [16] . It describes atomic repul- 
sion, due to the radiation pressure of scattered photons 
on nearby atoms, and can also be determined by a Pois- 
son type of equation 



V • Fr = n[r) 



(6) 



where gr is the atom scattering cross section. A detailed 
discussion of these forces and explicit expressions for the 
cross sections gr and gt can be for instance found in 
[HI US]. These expressions for the forces acting on the 
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where n and v are the density and velocity field of the 
gas respectively, Ft — Fmot + F c and P is the gas 
pressure. Here we have also defined the collective force 
F c = Fa + Fr , which is determined by the Poisson equa- 
tion resulting from equations ^ and ([6]), 



V-F c = Qn , Q = ((Jr- a L )(T L I /c 



(9) 



The system can then be regarded as a one component 
plasma where the electrostatic attractive force due to ions 
can be formaly replaced by the confining force, in such 
way that, for the unperturbed state n — no, we may 
represent the total force through a Laplace- like equation 



V • F T (n = n ) = 0. 



(10) 



In typical experimental conditions the repulsive forces 
largely dominate over the shadow effect, and the quantity 
Q is positive [41EIIH], which allows us to define the typical 
frequency 



U!p = 



Qnp 
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which is a straight-forward generalization of the well- 
known electron plasma frequency. Comparing this ex- 
pression with the usual definition of the plasma frequency 
Lu pe in an ionized medium, we conclude that neutral 
atoms behave as if they had an equivalent electric charge, 
as first noticed by [7], with the value q e ff = \ZtoQ, 
where eo is the vacuum electric pcrmitivity. The experi- 
mental value observed for this effective atomic charge is 
q e jj ~ 10 -4 to 10 -6 times the electron charge. In a typi- 
cal MOT experiment, we expect to have Uq w 10 10 cm~ 3 , 
M Pi 10~ 25 Kg, what yields a plasma frequency in the 
range of ujp/2n — q e ff /2ityJ m e jMu> pe f=a 100 Hz. It is 
clear from equation ( pT| that plasma like oscillations are 
only possible for Q > 0. Therefore, they cannot occur 
when the shadow force Q dominates over the repulsive 
scattering force ^ . 

In order to conclude the analogy between the fun- 
damental parameters between the plasma and the cold 
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atoms, one remark should be made concernig the motion 
of the center-of-mass 



R 



N 



flrrii/M, 



(12) 
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where M — X)iLi TO i- Therefore, for a typical neutral 
plasma with spherical geometry, it is a well know result 
that the center-of-mass oscillates at the so-called Mie fre- 
quency u>m — Wpe/V3) where the factor -\/3 arises from 
the spherical symmetry [TJ] and is an essential parame- 
ter in the description of resonances in clusters [T3] . How- 
ever, for the case of a spherical cloud of cold atoms, it is 
a simple task to verify that the center-of-mass obeys the 
following equation of motion 
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which states that the R oscillates exactly at the dipolc 
frequency lu^, corresponding to the generalized plasma 
frequency top in the unperturbed state n = no, as given 
by eq. (11). The reason for this difference remains in the 



fact that, in the case of the spherical plasmas, the restor- 
ing force depends on the charge balance between ions and 
electrons rii — n e , which on the other hand depends on 
the shape of the cloud. In the case of cold atoms, the po- 
tential is unequivocally determined by the laser and the 
magnetic field (say trapping) parameters, and therefore 
the shape and the atomic cloud plays no role. 

In what follows, we make use the set of fluid equations 
derived in this section for a quasi-continum medium and 
describe the main collective modes. We start from the 
infinite medium approximation, to emphasize the nature 
of the oscillations, and then introduce in Sec. IV the 
effect of the finite size of the cloud. 



III. PLASMA HYBRID WAVES 

We first assume oscillations that can be excited in the 
cold gas with a wavelength much smaller that its radius. 
The medium can therefore be assumed as infinite. We 
then assume that the equilibrium state of the gas is per- 
turbed by oscillations with frequency to and wavevector 
k. In the sense of linear response theory, we linearize the 
above fluid and Poisson equations, by defining perturba- 
tions around the equilibrium quantities 



n = no + n, F = Fq + 5F, v = vq + Sv. 



(14) 




FIG. 1: (color online) Normalized dispersion relation uo/up 
plotted against k\o- The red dashed line represents the 
acoustic assimptotic behavior of the hybrid waves. 



of simplicity. For the closure of the system of fluid equa- 
tions, one equation of state for the hydrodinamical pres- 
sure must be given. In this paper, we assume that P is 
given by an adiabatic equation of the form 



P(n) 



(15) 



where 7 represents the adiabatic constant. The impli- 
cations of this assumption will be stated in Sec. V, in the 
context of a kinetic description. Using the later together 
with equations ^ , ^ and ^ , we can easily obtain 
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where us can be identified with the sound speed 
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and Pq is the equilibrium gas pressure. 

We will now assume that the atomic equilibrium den- 
sity no is uniform, in consistence with the harmonic na- 
ture of the trapping potential, therefore neglecting the 
right hand side of equation ( 16 1. The influence of bound- 
ary conditions and inhomogencities on the collective os- 
cillations of the gas will be discussed later. Assuming 
a space-time dependence of the perturbations n and SF 
of the form exp(ifc • r — icot), with a complex frequency 
uj = uj r + iuji, we obtain for the dispersion relation and 
for the corresponding damping rate, the values 



„2„,2 



-a 



(18) 



Since the trapping force Fmot defines only the equilib- 
rium quantities and plays no role in the modes we are 

about to describe, we drop the subscript c for the per- dispersion relation reduces to J 1 — uj'p + k 2 u\, which is 

turbation in the collective force SF C in ( 14 1 for the sake formally identical to the dispersion relation of electron or 



In the limit of very small viscosity a -C Up, the later 
dispersion relation reduces to ui 2 = ojl 
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plasma waves in ionized media (also known as Langmuir 
waves) but where the electron thermal velocity Vthe = 
yfkBT e /m e , (T e and m e are the electron temperature 
and mass) is replaced by the sound velocity divided by a 
numerical factor us/V3- According to the experimental 
parameters refereed in the previous section, we estimate 
the sound velocity to value iig « y^kpT/M « 20 cm/s. 
This shows that the wave mode described by equation 
( [l8| contains elements of both electron plasma waves 
and acoustic waves. It possesses a lower cut-off, given 
by uj r — \/ up + 3a 2 /4, which is typical on an electron 
plasma wave, but its phase velocity tends to the sound 
velocity us and becomes weakly dispersive as an acous- 
tic wave. Its corresponding quasi-particles can therefore 
be seen as hybrid entities, somewhere between plasmons 
and phonons. The excitation of such modes in a typical 
MOT setup could be driven by modulating one of the 
six trapping laser beams. Such modulation should be 
nearly resonant with top, and hence one may excite hy- 
brid waves by sweeping the modulation frequency around 
100 Hz. Associated with the propagation of such waves, 
we expect to observe a periodic variation on the inten- 
sity of the luminosity. By measuring the period of this 
oscillation, one can identify the excited mode. 

The existence of such an hybrid mode is one of the 
main results of this paper. 



IV. MODIFIED TONKS-DATTNER 
RESONANCES 

The hybrid mode discussed above is only meaningful 
in infinite and homogeneous media. In physical terms, 
its dispersion relation can only be applied to waves that 
propagate locally, with wavelength scales much smaller 
than the inhomogeneity scale and the cloud dimensions. 
Let us consider now oscillations with a wavelength that 
is comparable with the size of the atomic cloud. In this 
case we can no longer neglect the boundary conditions. 
Going back to equation (16 1, we assume that the density 



except at the boundaries x — and x — L. Equations 



perturbations oscillate at a frequency u> as previously, but 
the corresponding spatial structure will be determined by 
the expressions 
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where the space dependent wavenumber k(r) is defined 

by 



(19 1 and (20 1 then reduce to a simple one-dimensional 



equation 
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Taking the boundary conditions n(0) = n{L) = 0, we 
obtain the following dispersion relation 



(22) 



where the quantum number v can take the values 
0,1, 2, 3,..., and the quantity = us /up is the Debye 
length for a cold neutral gas, in analogy with the plasma 
definition (where however the sound speed us is replaced 
by V^Vthe, as mentioned before). This defines a natural 
length above which plasma effects should be expected. 
Following the previously estimated values for ujp and 
us, we expect to observe a Debye length of the order of 
Xp, ~ 100 /im. In a typical MOT experiment, the radius 
of the cloud varies between a ~ 1 — 5 mm, yielding 
the relation Xp>/a <C 1. Hence, plasma-like effects are 
expected to take place in a cloud containing a moderate 
number of atoms. As a remark, we should stress out that 
there is no experimental evidence so far of any equation 
of state P{n), which may compromise the definition of 
the sound speed us- 



The relation (22) shows that the finite dimensions of 



the slab imply the existence of a series of resonant modes 
with an integer number of half- wavelengths. The cylin- 
drical geometry was considered, for the plasma case, in a 
famous paper by Parker, Nickel and Gould in 1964 |llj . 
but it is more natural here to consider a spherical geom- 
etry for the ultra-cold gas which, to our knowledge, was 
not derived for a plasma. We expect to find an infinite se- 
ries of resonances, similar to equation p2| , known as the 
Tonks-Dattner resonances [20j|2T]. For this purpose, we 
consider the internal oscillations in a spherical cloud with 
radius a in the homogeneous case, where Vn (r) = 0, for 
> r < a, for which analytical solutions can be found. 
These results remain qualitatively valid even for a more 
realistic density profile, as shown in Ref. [TT] , for the 
case of a inhomogeneous cylindrical plasma. The case 
of a non-uniform spherical plasma will be addressed in a 
future work. 

Performing a separation of variables, we can obtain 
solutions of the form 



fc 2 (r) = [c 2 -c4(r)]/4 (20) 

Before going into a more complex model, it is useful to 
consider the simple one-dimensional problem |19j . In the 
case of a uniform slab of cold gas, we have Vno = 0, 



h(r) = R{r)Y(9, </>) (23) 

where (r, 9, <f) are spherical coordinates. After separation 
of variables, we get the usual spherical harmonics for the 
angular part of the density perturbation 
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Y(6, <j>) = P 2 m (cos 0) exp(im<, 



(24) 



where P ; m (cos0) are the associated Legendre polynomi- 
als, I is a positive integer or zero, and |m| < Z. The radial 
equation resulting from (19 1 and (23 1 can be written as 
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By using a simple transformation of variables, x = kr, 
and S(kr) = y/kr R{f) this equation can be reduced to 
a Bessel equation 
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The solutions with regular behavior at the origin x = 
are therefore given by Bessel functions of the first kind, 
•h+1/2 {%)■ From this we conclude that the Tonks-Dattner 
modes in a spherical homogeneous cold atom cloud are 
determined by 



h(r, t) 



V hi{t) Jl+1 '^ Pl n {cos 6) exp(im^ 
, vkr 



(27) 
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FIG. 2: (color online) Normalized modes u^/up plotted 
against the nodes r^j /a of the radial solution for the density 
perturbation h, for Au/a = 0.1 and 1 < v < 10. Blue squares 
(1 = 0), red circles (Z = 1), violet stars (1 = 2), yellow triangles 
(1 = 3), green open circles (1 = 4). The full line is plotted at 
I = and scales as l/u 2 . 



where fii(t) have small amplitudes such that \fii\ <C no. 
The mode frequencies can be obtained by remarking that 
n should vanish at the border r — a. This implies that the 
allowed values for k have to obey the condition k — z Ut i /a, 
where z„ t i represents the fth zero of the Bessel function of 
order (I + 1/2). We are then led to the mode frequencies 



UJp 



1 + Z V A 



(28) 



Comparing with the rectangular case of equation ( 22 ) 



we see that the allowed eigen-frequencies for a spheri- 
cal cloud now depend on two quantum numbers v and 
I. But, in contrast with the similar quantum mechani- 
cal solutions for hydrogen like atoms, we have no hier- 
archical relation between these quantum numbers. The 
normalized radial profiles for the lowest order solutions 
are illustrated in figure ([2]), which is formally similar to 
the spectrum that was presented in the past by Dattner 
for the case of a plasma cylinder |21J . Also, a recent ex- 
perimental on ultracold plasma was published showing 
the evidence of Tonks-Dattner modes excited by a radio- 
frequency electric field [23] , 

In a more realistic description, the present rigid 
(Dirichlct) boundaries will eventually have to be replaced 
by soft boundaries and a generic density profiles no(r) 
must be assumed, for which numerical solutions have to 
be found E21 ■ 



V. NONLINEAR OSCILLATIONS DRIVEN BY 
DIPOLAR RESONANCES 

Due to the intrinsic nonlincarity of the fluid description 
the collective behavior of the medium, it is possible to 
couple dipole oscillations of the cloud with the hybrid 
modes. Going back to the fluid equations, and assuming 
an oscillating mean velocity of the form vq sin(w c ;t + <f>) 
for the center of mass, we obtain 



-Wo 
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dn 
] l)t 
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(29) 

Using the factorization h(r,t) = A(t)h(r), the coupling 
between the center of mass and the hybrid modes can be 
approximately described by the canonical Mathieu equa- 
tion 



d 2 A 

W 2 + 

where we have used 



[A + 2ecos(2r)]i = 0, 



(30) 
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It is well known that such an equation has unstable re- 
gions. For |e| <C 1, the first transition to the instability 
zone occurs for A ~ 1+e. We therefore expect to observe 
an instability of the hybrid and Tonks-Dattner modes, 
driven by dipolc oscillations for 
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vq ■ V In h 



(32) 



This simple discussion demonstrates the existence of very 
interesting nonlinear collective phenomena in a cold atom 
gas. But in the present work we are mainly focused on 
the linear properties of the medium. A more rigorous and 
detailed study of nonlinear coupling between dipole and 
Tonks-Dattner oscillations will be left to future work. 



VI. KINETIC DISPERSION RELATIONS 

The analysis of collective oscillations in a cold gas can 
be refined by using a kinetic description based on the 
classical probability distribution function W{f, q, t). This 
description will allow us to include resonant kinetic pro- 
cesses, which enhance the energy exchanges between part 
of the atomic population and the hybrid modes excited 
in the medium. Moreover, a kinetic approach avoids the 
use of an equation of state for the hydrodinamical pres- 
sure P. It is rather determined by the statistics, as we 
will see in the following discussion. In the presence of 
diffusion, the probability distribution function W(r,q,t) 
obeys the Fokker-Planck equation [21] 
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(33) 

where the total force F to t includes the radiative and the 
damping forces, and the diffusive tensor is due to 
the fluctuations of the radiative force and spontaneous 
emission, as discussed by several authors [24, 25 . Here, 
for simplicity, and because we want to focus on the os- 
cillating modes, we neglect the diffusion term. Diffusion 
effects will be not completely ignored, rather they will 
reappear later in a different context. We are then led to 
a kinetic equation of the Vlasov type, with a damping 
correction, as given by 




v- V 




W = -a(v)W (34) 



where v — fiq/M is the atom velocity, and the collective 
(shadow minus repulsive) force F c is determined by a new 
Poisson equation, which can be written as 



(35) 



V -F C = Q / W(f,v,t)dv 



This equation is obvious identical to H , because the inte- 
gral is nothing but the density n(r, t), the normalization 
of the Wigner function. On the other hand, by taking 
the average of the first two momenta of the kinetic equa- 
tion ^34]), we will be able to derive fluid equations such 
as (|9|-( 10 1. Notice that, the parameter a appearing in 



the fluid equations in an averaged value of the quantity 
a(v) appearing in the wave kinetic equation (34 1. In or- 



der to focus our attention on the purely kinetic processes, 
we will assume that a = but, at the end of this sec- 
tion, we will discuss the the influence of a finite value for 
this parameter. We now consider some equilibrium state 
Wq(v) and assume a sinusoidal perturbation, such that 
SF and W evolve in space and time as exp(ifc ■ f — iujt). 
In what follows, we drop the subscript c in the collective 
force, in consistence with the previous calculations. After 
linearization, the two previous equations reduce to 



W= - 



i 5F-dW a /dv 
(uj — k ■ v) 
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ik-SF = 



W(v) dv 



(36) 

From here we get the dispersion relation for collective 
cold atom oscillations with frequency u and wavevector 
fc 



Q r k- dWo/dv 
~Mk* J {oj-k-v) 



(37) 



This is similar to that of electrostatic waves in unmagne- 
tized plasmas, and can be rewritten as 1 + x(w, fc) = 0, 
where the quantity x( w j k) is the susceptibility. In order 
to understand the physical implications of such a disper- 
sion relation, let us consider first a simple mono-kinetic 
atomic distribution, of the form Wq[v) = n^8(v — v ), 
corresponding to a beam of atoms with density no and 
velocity vq. In this case, equation d36| reduces to 



Qn 

M(LO-k-V ) 2 



= 



(38) 



This is nothing but the Doppler shifted plasma oscil- 
lations discussed above. For vq = 0, this reduces to 
uj = ujp = (Qtiq/AI) 1 / 2 . Second, we consider the case of 
two atomic populations, one at rest with density n , and 
another moving with velocity v\ with density ri\ < no. 
This is described by Wo(v) = no5(v) + ni5(v — V\). The 
corresponding dispersion relation is given by 



Qn 



MoS 1 M(ui - k ■ vif 



(39) 



The interest of this new dispersion relation is that it leads 
to an instability, which means that the cold atoms can 
oscillating spontaneously at a frequency close to uj p. The 
energy source driving such oscillations is provided by the 
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atomic beam with density n\ . This can be seen by assum- 
ing a complex frequency to — uj r + iu>i. The maximum 



growth rate predicted by equation (|39| corresponds to 
the resonant case where u) r 

by 



= k 



Vi, 



and is determined 



Ulp 



2n 



1/3 



(40) 



This is very similar to an electrostatic beam-plasma in- 
stability, but can also be seen as the dynamical analogue 
of the collective atom recoil laser (carl) [26], where the 
incident laser beam is replaced by the atomic beam, and 
where the resulting stimulated emission is not made of 
photons, but of hybrid quasi-particles. The critical re- 
gion, where unstable modes may show up, is given by 



K < (1 + iV 2 / 3 ) 3 / 2 , (41) 

where we define the dimensionless quantities K = k ■ 
vi/cjp, Q = lu/lu p and N = rii/riQ. In Fig 3j we plot the 
growth rate T(K) = Ui/u) p for the unstable region given 
above and the real part of the spectrum fl r = u) r /u) p , 
both given by the roots of the dispersion relation [39] We 
observe that the stable modes bifurcate exactly at the 
vanishing point of the unstable ones. 

Finally, we consider a generic equilibrium quasi- 
distribution Wq(v). The atomic susceptibility can be 
split into its real and imaginary parts, \ = Xr + 
By using the Kramers-Kroenig relation between x.r and 
Xi, and assuming that most of the atoms have velocity 
smaller than the phase velocity uj/k, which is a very plau- 
sible assumption for a gas of ultra-cold atoms, we obtain 




FIG. 3: Normalized roots of Eq . J39| in terms of the dimension- 
less quantities fl = lj/uj p , K = k ■ vi/ui p and N = ni/no. (a) 
Dimensionless growth rate T(K) = +£li(K) (unstable solu- 
tions) and (b) dimensionless modes Q r (K) (stable solutions). 
In both plots, N = 1 (full black line), N = 1/2 (dashed violet 
line) and N = 1/3 (dotted blue line) . 



Xi(w,fc) = ™ 



Q fdG 



Mk 2 V dv 



j/k 



(42) 



where v is the parallel component of the atom velocity, 
and we have identified the sound speed with the integral 



This means that the hybrid modes are essentially a one- 
dimensional process, if we remind that 7 = (2 + d)/d, 
which is a well known result for the electron plasma 
waves. This justifies a posteriori the hybrid character 
of this oscillations. Back to eq. (39 1, we can easily ob- 
tain the dispersion relation, by using 1 + Xr(^r,k) = 0, 
which coincides with equation (18), and the wave damp- 
ing defined by the expression 



4 = — / W (v)v 2 dv= — 



G {v)v 2 dv (43) 



The quantity Gq(v) introduced here is the average of 
the Wigner distribution over the perpendicular velocities. 
The later result avoids the postulation of the equation of 
state (15 1. However, since the hydrodynamical pressure 



P(r, v, t) = -Mn(r, v, t) ((v 2 ) - (v) 2 ) (44) 
3 



is a statistical variable we can make use of eqs. ( |43| and 
( 17 1 to write 



7 



(45) 



(dxr/duj)^ 



7T Q f dG 
u> Mk 2 \ dv 



oj/k 



(46) 



This is a non dissipative wave damping, which is not re- 
lated with any increase of the entropy of the physical 
system, which is the so-called Landau damping. It de- 
scribes the resonant interactions between the wave and 
the atomic population with has a parallel velocity nearly 
equal to the wave phase velocity. Usually, for a thermal 
equilibrium distribution W$(v) this quantity is negative, 
and corresponds to wave damping. But the sign of u>i can 
change for a non-thermal distribution, eventually leading 
to wave instability and wave growth. 
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We can take a step further in the kinetic description 
of the collective oscillations in the cold atom cloud, and 
consider a broad spectrum of fluctuations, described by 
the total wave intensity 



J(i)= / i(M) 



dk 
'(27r)E 



(47) 



where the spectral intensity is defined by I(k, t) — 
W* (k, t)W(k, t). Following the usual steps of the plasma 
quasi-linear theory, and adapting it to the present con- 
text, we can say that, each spectral component behaves 
in accordance with the above description, and evolves in 
time according to the equation 



dt 



I(k, t) = 2co,(k, t)I(k, t) + S(k, t) 



(48) 



where S(k,t) is any given source term, and the total 
damping rate 7& (t) slowly evolves in time due to the slow 
time evolution of the equilibrium (or quasi-equilibrium) 
distribution Wo(v,t), which can only be considered con- 
stant in a short time scale. The temporal evolution of 
Wo(v,i) under the influence of the fluctuation spectrum 
is determined by a diffusion equation of the form 



where the new diffusion tensor D, associated with the 
collective oscillations, is determined by 



D(tf,t) 



Ldp 



I(k,t) 



k®k 



dk 



(w-jfc-jfc) ( 2 ^) 3 



(50) 



Comparing this with our previous kinetic equation ( 33 ) 



it can be seen that the existence of a collective spectrum 
of oscillations introduces an additional diffusion effect in 
atomic velocity space, which tends to prevent the atomic 
cooling process. As we have noticed, these results are 
only valid in the limit of a negligible viscosity parameter, 
a — ► 0. A finite value of a will have two distinct con- 
sequences. First, it will lead to the damping coefficient 
already stated by equation (18 1, adding to the purely ki- 
netic Landau damping. Second, it will broaden of the 
Landau resonance appearing in the dispersion relation 
( 37 1 , and in the diffusion coefficient ( 50 1 therefore reduc- 



ing the efficiency of the resonant atom-collective wave in- 
teractions associated with the Landau resonance and to 
the quasi-linear diffusion. It will therefore compete with 
the kinetic effects described in this section. The com- 
bined influence of viscous and kinetic damping is outside 
the scope of the present work and would deserve a sepa- 
rate investigation. 



VII. CONCLUSIONS 

In this work we have used both fluid and kinetic equa- 
tions to describe the collective oscillations in a cloud of 
neutral atoms confined in a magneto-optical trap. Our 
approach is based on a simple but physically relevant 
model for the forces acting on the cold atoms, which has 
been well verified by the experiments, and can be de- 
scribed by a Poisson equation similar to that describing 
electrostatic interactions [H [SJ [TS]. Once the physical 
picture was established, we started by setting the basic 
equations and parameters. We have also shown that the 
existence of internal forces leads to the existence of col- 
lective waves which have an hybrid character, with prop- 
erties that are common to both electron plasma waves 
and acoustic waves. The similarities and differences of 
the cold atom gas oscillations and those of a plasma were 
discussed. 

Taking into account the finite size effect in a cloud of 
cold atoms, we have shown that internal resonances of 
these hybrid oscillations inside the cloud can be excited. 
These modes should be called Tonks-Dattner resonances, 
in analogy with similar plasma physics effects, since they 
are shown here also to exist in a neutral magneto-optical 
trapped gas as well. Previous models developed for plas- 
mas were limited to planar and cylindrical geometries, 
and they were extended to the spherical geometry which 
is more relevant to the cold atom clouds. We have ex- 
tended our discussion of these oscillations to the non- 
linear regime, where coupling between dipole oscillations 
and hybrid waves can take place. Our approximate de- 
scription suggests that such a nonlinear coupling can lead 
to the destabilization of hybrid waves or Tonks-Dattner 
resonances, driven by dipolar oscillations of the cloud 
center of mass. 

Our analysis proceeded with a kinetic approach. Us- 
ing this more refined description of oscillations, we were 
able to derive more general dispersion relations where 
non-dissipative Landau damping was included. We have 
limited our discussion to hybrid wave modes with wave- 
lengths much smaller than the typical dimensions of the 
atomic cloud. But, the results can easily be extended 
to the confined Tonks-Dattner resonances. Finally, we 
have established a quasi-linear kinetic equation, show- 
ing the occurrence of diffusion in atomic velocity space. 
This diffusion effect is a direct consequence of the col- 
lective fluctuation spectrum, and implies the existence 
of new collective processes preventing the occurrence of 
laser cooling. 

In the present work we have explored the similarities 
of the could atom cloud with a plasma, which can be as- 
sociated with the existence of an effective electric charge 
for the neutral atoms. The resulting wave modes however 
are not identical to the plasma wave modes, but show an 
hybrid character. We hope that this work will motive 
future experimental and theoretical work on the collec- 
tive oscillations in cold atom traps, and will contribute 
to launch cold atom research in new directions. 
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